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Abstract 

Although marginally more complicated than the traditional Laplace sum-rules, Gaussian sum-rules have the 
advantage of being able to probe excited and ground states with similar sensitivity. Gaussian sum-rule analy- 
sis techniques are applied to the problematic scalar glueball channel to determine masses, widths and relative 
resonance strengths of low-lying scalar glueball states contributing to the hadronic spectral function. A feature 
of our analysis is the inclusion of instanton contributions to the scalar gluonic correlation function. Compared 
with the next-to-leading Gaussian sum-rule, the analysis of the lowest-weighted sum-rule (which contains a large 
scale-independent contribution from the low energy theorem) is shown to be unreliable because of instability under 
QCD uncertainties. However, the presence of instanton effects leads to approximately consistent mass scales in the 
lowest weighted and next-lowest weighted sum-rules. The analysis of the next-to-leading sum-rule demonstrates 
that a single narrow resonance model does not provide an adequate description of the hadronic spectral function. 
Consequently, we consider a wide variety of phenomenological models which distribute resonance strength over a 
broad region — some of which lead to excellent agreement between the theoretical prediction and phenomenological 
models. Including QCD uncertainties, our results indicate that the hadronic contributions to the spectral function 
stem from a pair of resonances with masses in the range 0.8-1.6 GeV, with the lighter of the two potentially having 
a large width. 



1 Introduction 

Mass predictions for scalar (0++) glueballs extracted from QCD sum-rules have been problematic mainly due to 
discrepancies between analyses which are sensitive to the low-energy theorem for gluonic correlation functions and 
those which are insensitive to this quantity H, ||, ^ ^, ^ . Such a discrepancy would be indicative of two widely- 
separated states, a result which has already been seen to occur in explicit two-resonance analyses of Laplace sum-rules 
even in the absence of mixing with quark scalar resonances 1^, |^. However, there exists substantial evidence that these 
discrepancies are resolved by the inclusion of instanton \f\ effects in the Laplace sum-rules for scalar glueballs [B[ |^] . 

Recently, techniques for using Gaussian sum- rules Q to predict hadronic properties have been developed f^. In 
particular, these methods concentrate on normalized Gaussian sum-rules which are independent of the finite-energy 
sum-rule constraint which is central to the original heat-evolution studies of the Gaussian sum-rules. Advantages 
of this approach compared with Laplace sum-rules include enhanced sensitivity to the hadronic spectral function 
over a wide range of energy. In this paper, these techniques are employed and generalized in an effort to obtain 
resonance parameter predictions for scalar gluonium. Furthermore, the formulation of these sum-rules is extended 
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to include Gaussian kernels weighted by integer powers, 



-p{t)dt , fc > -1 (1) 

TT 
to 

where p{t) is a hadronic spectral function with physical threshold to- Similar to the Laplace sum- rules, the low- 
energy theorem (LET) (see (^ below) for scalar gluonic currents enters only the fc = — 1 Gaussian sum-rule, 
and instanton contributions to the correlation function serve to mitigate the discrepancy between the k = —1 and 
fc > — 1 sum-rules. However, theoretical uncertainties associated with the instanton and LET parameters are shown 
to be overwhelming in the fc = — 1 sum-rule, rendering it unsuitable for phenomcnological analysis. Thus the fc = 
Gaussian sum-rule is the focus of our detailed predictions for scalar glueballs, including an estimate of theoretical 
uncertainties. 

We show that the Gaussian sum-rules of scalar gluonic currents contain signatures that the hadronic spectral func- 
tion is distributed over a broad energy range, (0.8-1.6) GeV including QCD uncertainties, and that this distribution 
most likely consists of two separate resonances. Since these sum-rules probe the gluonic content of hadronic states, 
they are sensitive to the glueball component of the observed scalar mesons which in general could be glueball-quark 
meson mixtures. Thus our results are relevant to the interpretation of the scalar isoscalar resonances in this region: 
the /o(400-1200), /o(980), /o(1370), /o(1500), and /o(1710) 

In the next section, Gaussian sum-rules for scalar gluonic currents are developed. In Section]^, we develop analysis 
techniques and employ them to analyze the Gaussian sum-rules using a variety of phenomcnological models. The 
results of this phenomcnological analysis, including theoretical uncertainties, are consolidated in Section ^, and a 
summary of our results is contained in Section 0. 
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2 Scalar Glueball Gaussian Sum-Rules 

The most important quantity in any sum-rules approach to determining hadron properties is the correlation function 
for the particular channel under inspection: 

n(Q2) /"d4a;e'9-"(rj|T{J(x), J(0)}|r!) , = -q^ (2) 



where \Vt) is the QCD vacuum state, T is the time-ordering operator, and J{x) is that current which corresponds to 
the quantum numbers of interest. In this paper, we wish to focus on scalar glueballs and so we choose the following 
current: 

J = -^/3(a)GVG'^^'^ (3) 

which is renormalization-group invariant in the chiral limit of nf massless quarks. The gluon field strength tensor 
G^^j, is defined by 

= d^Al - d.AI + gr'^A^Al (4) 
and P{a) is the QCD beta function describing the momentum scale dependence of the strong coupling parameter a 
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^o-^-Q^f - = , ... . (6) 

From the asymptotic form and assumed analytic properties of (^) follows a dispersion relation with three sub- 
traction constants 



n(g^)-n(o)-Q^n'(o)--g^n"(o) = -:^y^ ^a^^l^.^ dt , q'>o (7) 
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where p{t) is the hadronic spectral functiorj^ with physical threshold Iq. The spectral function p{t) is related to 
a physical process and is thus determined phenomenologically. In contrast, ll{Q^) is calculated theoretically from 
QCD, and the constant 11(0) follows from the low-energy theorem 

n(0)^ lim U{Q')^^(J) . (8) 

For these reasons, we shall refer to the left-hand side of as the theoretical side and the right-hand side as the 
phenomenological side. In this regard, eqn. (|^) serves to relate theory to phenomenology, and, in principle, could be 
used to predict the properties of hadrons from QCD. 

However, as it stands, eqn. (|^) is not actually that well-suited to this task. For instance, although the constant 
n(0) is determined by the low-energy theorem (|8|), the constants n'(0) and n"(0) are not. Further, the theoretical 
calculation of n(Q^) contains a field theoretical divergence proportional to Q"*. In addition, from a phenomenological 
perspective, the integral on the right-hand side of (Q) is far too sensitive to the high energy behaviour of p{t) to 
effectively probe low-lying resonances. 

To circumvent these shortcomings, we consider the one-parameter family of Gaussian sum-ruleg|^ 

G.js, r).^B[^'^ - - - ^^^''"(-^ + } , k>-l (9) 

with the Borel transform B defined by 



A^/N=At 

Applying definition (||) to both sides of (0) alleviates the difficulties surrounding (Q): the infinite number of deriva- 
tives in (|l0| ) annihilate the unwanted low-energy constants and the field theoretical divergence contained in Il{Q^). 
Furthermore, as we shall see, a key feature of the resulting sum-rules is the introduction of a Gaussian weight 
factor to the integrand on the phenomenological side of @). This serves to suppress contributions from p{t) away 
from the Gaussian peak — a desirable situation considering that we wish to extract information concerning low-lying 
resonances. 

Let us first consider (^) as applied to the theoretical side of (0). As noted previously, the low-energy constants 
n'(0) and n"(0) are annihilated by the Borel transform; however, the constant 11(0) does produce a contribution 
unique to the case k = —I. Using the following identity: 



B 



A2 



= — (-a)"exp| — ) for?i>0 , (11) 
4r V 4t ' 



it is trivial to show that the contribution to the theoretical side of the Gaussian sum-rules devolving from the 
low-energy constant is given by 



To proceed further, however, we must settle on a specific form for the scalar glueball correlator Il{Q^). We choose 
to partition the correlator into the following sum of qualitatively distinct terms: 

jjQCD^g2) = np°''*(Q2) + n™"'i(g2) + n'"^'(Q2) , (i3) 

where the superscript QCD signifies that (|l^) is a theoretical approximation to the true correlator. The first two terms 
in (^3|) devolve from the operator product expansion of the current The quantity lU"^'^*^ {Q'^) is the contribution 

'^In the literature, p(t) is often denoted by Imn(i). 

^This definition is a natural generalization of that given in [pi. To recover the original Gaussian sum-rule, we simply let fc = in (P). 
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from ordinary perturbation theory whereas 11'^°"'^ (Q^) is the result of nonzero vacuum expectation values of local 
gluonic operators (condensates). For three colours and three flavours of massless quarks {rif — 3), IU'^'^*'{Q^) is given 
at three-loop order by |T2] 
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02 = -10.1250 



(14) 

(15) 



where we have omitted the aforementioned field theoretical divergence as this term is annihilated by the Borel 
operator (|l0|). Incorporating into 11'^°"'^ (Q^) next-to-leading order [|l^ contributions^ from the dimension four gluon 
condensate (J) and leading order Q contributions from gluonic condensates of dimension six and eight 



(Ce) — (gfabcG'^^GlpG'p^) 

(Os) = 14((a/,b,G;,G^,)') - {{afabcG'J.^Gl^) 



yields 



jjcond^g2) = 



&0 = in— 

TT 
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(16) 
(17) 

(18) 
(19) 



The final term on the right-hand side of (|13|) is a contribution arising from direct instanton effects. Note that 
in decoupling this term from perturbation theory and the condensate contributions, we have tacitly assumed that 
interference between classical and quantum fields is small. Further, we also assume that the dominant contribution 
to n'"'^'^(Q^) comes from BPST single instanton and anti-instanton solutions |^] and that multi-instanton effects are 
negligible iQ. With these provisions, we have ||l|, H, [l5|| 



K2 (pVO^)]" dnip) 



(20) 



where K2{x) is the modified Bessel function of the second kind of order two (c./. ||lq| ), p is the instanton radius, and 
n{p) is the instanton density function. 

Before substituting (13) into (^), it is convenient to first simplify (^ by employing a particularly useful identity 
relating the Borel transform (|l^ to the inverse Laplace transform jsj 



(21) 



where, in our notation, 



1 l-a+ica /A2\ 



(22) 



with a chosen such that all singularities of / lie to the left of a in the complex A^-plane. (Note, there is no loss of 
generality in assuming that a > 0.) Rewriting (^) using (^) gives 
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(23) 



^The calculation of next-to-leading contributions in |l3| have been extended non-trivially to rtj^ = 3 from rif = 0, and the operate 
basis has been changed from (^aG^) to (J). 
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If in the first integral above, we make the substitution w = ~s — iA and in the second, we make the substitution 
w = — s + iA, then (p3|) reduces to 



Gfc(s,T) 



1 I 

Ay/TTT 2711 
I 1 
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{-w)^ exp 
[~w)^ exp 



-(s + 
47 



n(w)dw + 2 / (-w)''exp 
n(i(;)dw 



47 



n(w)dw 



where Fi and r2 are two parabolas (depicted in Figure y) in the complex w-plane defined by 

Fi = -s-i(a2+x2)i/4exp 
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i(a2+x2)i/4 
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(24) 
(25) 

(26) 
(27) 



for all X e R. 

Now, we must substitute ( |l3|) into (^5|) and calculate the resulting complicated integral. Towards this end, it is 
advantageous to consider the closed contour C{R) depicted in Figure |[ Our expression for the correlator (|l^) is 
analytic in the complex w = Q^-plane except for a branch cut along the negative real semi-axis originating from a 
branch point located at the origin. Consequently, 
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(28) 
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where Fi(i?) and F2(i?) are respectively those portions of the contours Fi and F2 (see ( |26| ) and ([27])) lying in the 
interior of a circle of radius R centered at — s. For large i?, the integral over F3 + F4 + F5 approaches zero and the 
contours f i(i?) and f 2(-R) approach Fi and F2. Therefore, by rearranging (p9[), recalling (p5[), and taking appropriate 
limits, we get 
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where 



ImnQ^°(rt = lim 



2i 



(31) 



(32) 



Eqn. (^) is our final expression for the contribution to the fc-th Gaussian sum-rule of scalar gluonium stemming 
from the correlator ([T^). Later in this section, however, we do evaluate ( ^l|) for the specific cases k = —1,0 (see (^) 
and(g|)). 

We must now consider the phenomenological contribution to the Gaussian sum-rules. Substituting the right-hand 
side of (0) into (||) and again making use of the identity ([III), it is simple to show that 
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In sum-rules analyses, it is customary to approximate the spectral function p{t) using a "resonance(s) + continuum" 
ansatz. In this model, hadronic physics is (locally) dual to QCD above the continuum threshold sq, and so we write 



p{t) = phad(^) +9{t-so) ImnQCD(t) 
where 9(t) is the Heaviside step function. (We shall have much more to say concerning p^'^'^(t) in Section ||.) 



(34) 



Substituting 
contribution 



Gr'(s,r,so) 



into (|33| ) and comparing the result to the theoretical expression ( plf) shows us that the continuum 

{s-tr 
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is common to both; therefore, we define 



G^^°(s, T, so) ^ GQ°°(s, r) - Gr\s, r, so) 
Gr^(J,T) ^ Gf <="(!, r) - Gr\s,T,so) 



and write (recall the low-energy contribution (|T^ unique to the k — —1 sum-rule) 



GQ?°(J,r,so)-f 
GQ^°(J,r,.o) = 
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(35) 



(36) 
(37) 



(38) 
(39) 

(40) 



We note that (38-40|) have meaning for s < since this simply represents a Gaussian kernel whose peak lies outside 
the t > Q physical region, and hence only the Gaussian tail extends into the physical region. As will later be shown, 
the QCD results scale with r through the renormalization-group equation, and hence QCD also presents no obstacles 
to considering s < 0. Indeed, the seminal work on Gaussian sum-rules considered symmetric and antisymmetric 
combinations (s, t) — G (s, r) ± G (— s, r) which implicitly employs Gaussian sum-rules with negative s. 

In this paper, we focus exclusively on the fc = 0,-1 Gaussian sum-rules. Substituting ( p^ into ( |36| ) and 
recalling (^ij) givesQ (for details on simplifying the relevant integrals in (^), see [|[ ^) 
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The given result is valid to leading order in the condensates. 
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where J2(x) and Y2{x) are Bessel functions of order two of the first and second kind respectively (c./. [p^). 
Renormalization-group improvements of and amount to replacing the strong coupling constant 

a (contained in the coefficients ( p^ ) an d (|f9| )) by the running coupling a{i>'^) at the renormalization scale — ^/t. 
At three loop order with rif = 3 in the MS renormalization scheme, we have fv\\ 

a[u^) _ 1 ^ilogL , 1 .-2 2 



[^2(log2L-logL-f)+^2] (43) 



with Ajgg- « 300 MeV for three active flavours, consistent with current estimates of a{M^) ^ and matching 
conditions through the charm threshold . 

The normalization of the Gaussian sum-rules is related to the finite-energy sum-rules (FESRs) |20| as can be seen 
by integrating (|3|) and (^9|) with respect to s to obtain 



j G'^'f{s, r, so)dJ + n(0) = i/^^(Odi (45) 

— OO Iq 
OO OO 

J G^™{s,T,SQ)ds^^ J t''p^'^'^{t)dt , A:>0. (46) 

— OO to 

We recognize the quantities on the right-hand sides of ( ^5| ) and ( ^6| ) from the definition of the (FESRs) 

OO 

Fkiso) = l J dtt'^p'^^it) , (47) 

to 

where Fk represents a QCD prediction. Thus we see that the overall normalization of G^*""^ (or G?^'^ -I- n(0)) is 
constrained by the finite-energy sum-rules. 

This result is not surprising in light of the seminal work on Gaussian sum-rules which established the significance of 
the FESR constraint by considering the evolution of the Gaussian sum-rules through the diffusion equation ||] . ft was 
found that this "heat-evolution" of the resonance plus continuum model would only reproduce the QCD prediction 
in the asymptotic regime if the continuum sq was constrained by the lowest FESR. Hence, the normalization of 
the Gaussian sum-rules, which is constrained by the FESR, should be removed by defining normalized (unit-area) 
Gaussian sum-rules 



TV?™ (s, T, so) = acT, (48) 

{s,t,sq)= , k>0 (49) 

where the n-th moment of Gk is given by 

OO 

Mfe,„(T,so)= / s"Gfc(s,r,so)ds , n = 0,l,2,... . (50) 



^The role of higher-loon_effects on the extraction of a (Mr) and the subsequent impact on A-j^ has been investigated through Pade 
approximation techniques llSl. 
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Note that, for the sake of notational convenience in subsequent sections, we have absorbed the low-energy theorem 
contribution into the definition of 7V*^{^'^ (sec ( p8| ) above). This ahows us to write our final version of the normalized 
Gaussian sum-rules of scalar gluonium as 



ivQ^°(s,r,.o) = 



At 



1 fOO , 



k> -I 



(51) 



Before describing analysis methods for the Gaussian sum-rules, we specify the numerical values for the QCD 
parameters appearing in (|l3|) that will be employed in our analysis. For the dimension four gluon condensate, we 
make the assumption that (J) w [aG"" ^^,0°-^^) and then employ the most recently updated value pl| 

(aG'^^^G"'''') (0.07 ±0.01) GeV^ . (52) 

The dimension six gluon condensate can be related to (aG^^^G"'^'') using instanton techniques (see |l|, [2^ ) 

(Oe) = (0.27 GeV2)(aG%G'''''') . (53) 

Further, by invoking vacuum saturation in conjunction with the heavy quark expansion, the authors of [ p3| have 
related the dimension eight gluon condensate to {aG"" ^^G""^^) through 



(54) 



Regarding the instanton contributions, we shall employ Shuryak's dilute instanton liquid model |p4| in which 

n{p) ^ ncS{p - pc) (55) 

with 



1 



= 8.0 X lO"'' GeV'' and p^ = — GeV^ 



0.6 



(56) 



3 Analysis of the Scalar Glueball Gaussian Sum-Rules 

In most sum-rules analyses, it is necessary to make some assumptions concerning the hadronic content of the spectral 
function (|3^). A common choice in the literature is that of a single narrow resonance model. However, the wealth 
of scalar states with masses below 2 GeV |l^, some of which are quite broad, certainly raises a question as to the 
validity of such an assumption, and suggest that a more general model allowing for a distribution of the resonance 
strength would be more suitable. In the following subsections, we show that this is indeed the case. Using appro- 
priate generalizations of the analysis techniques developed in , we demonstrate that the single narrow resonance 
model provides an inadequate description of the hadronic content of the scalar glueball correlator, whereas certain 
distributed strength models lead to outstanding agreement between the theory and phenomenology parameterized 
by i^. 

3.1 Single Narrow Resonance Model 

In the single narrow resonance model, p^^'^{t) takes the form 



(57) 



where m and / are respectively the resonance mass and coupling. With such an ansatz, the normalized Gaussian 
sum-rule ( ^T| ) becomes 



7VQC°(s,T,5o) 



2\2 - 



: exp 



(s — m ) 
47 



(58) 
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The phenomenological side of ( pq) a dmits an absolute maximum (peak) located at s = m^, independent of 
r; therefore, the theoretical side o05q) should mimic this behaviour provided that the single narrow resonance 
model ( p7| ) is actually an adequate description of hadronic physics below the continuum threshold {i.e. heavier states 
are weakly coupled enough to be absorbed into the contimmm). Defining Spcak(''", sq) by the condition 



_^^QCD 

ds ^ 



iV^^^°(spoak(^:So),r,so) =0 , (59) 



and denoting by {t„}^^q an equally spaced partition of the r interval of interest [rj, rf] (we elaborate on this interval 
shortly), we define the following x^-function 

N 

X^(S0, m^) = ^ [speak(T„, Sq) - TO^] (60) 
n=0 

as a measure of the difference between the theoretical peak position and the phenomenological peak position. Mini- 
mization of ( |60| ) with respect to Sq and then provides us with values for these two parameters which correspond 
to the best possible fit between the theory and phenomenology as represented through (^8|). Lastly, the optimizing 
condition 

-^X'(so,m2) = (61) 
allows us to write as a function of so whereby reducing ( |60| ) to a one-dimensional minimization problem: 

N 

X^{sq) = ^ [spoak(T"„, So) - in^iso)] (62) 

n=0 

with 

N 

rr? 

Thus, in a single narrow resonance analysis, we first minimize ( |62| ) with respect to sq to determine an optimum choice 
for the continuum threshold parameter and then substitute this value into ( |6^ ) to obtain the best fit resonance mass. 

Upon obtaining these optimized parameters, there exists criteria that can be used to assess the validity of 
this phenomenological model. For instance, if the single resonance analysis is a reasonable approach, plots of the 
theoretical Gaussian sum-rules (the left-hand side of (^8|)) and plots of the phenomenological Gaussian sum- rules 
(the right-hand side of (^)) should coincide (to a large degree). Significant deviation of one from the other may be 
indicative of an inadequate phenomenological model. On a more quantitative note, consider the following combination 
of moments (pOh (where we suppress the explicit dependence on r and sg)- 



(so) = -^fl^ 5Z ■5poak(T„, So) . (63) 



,2 _ 



(64) 



In the single narrow resonance model we should find 

al^2T^Q , (65) 

and hence a substantial deviation of cr^ — 2r from zero signals a failure of the single narrow resonance model. 

To proceed with the analysis, however, we must first choose our region of interest [rj, Tf] needed in the definition 
of the function (|60|). There are a number of factors to be considered in selecting this interval. The lower bound 
Tj must be large enough such that the condensate contributions do not dominate perturbation theory and also such 
that the leading omitted perturbative term in the expansion for the running coupling (^ ) is small. Therefore, in 
accordance with these criteria, we choose a lower bound of Tj > 2 GeV^. To choose an appropriate upper bound on r, 
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we first note that the Gaussian Icernel has a resolution of v2t. It is important to the analysis that the Gaussian sum- 
rules employed have a resolution less than the non-perturbative (hadronic physics) energy scale involved: roughly 
2-3 GeV^. This fact motivates an upper bound of Tf < 4 GeV^. Therefore, in this and in all subsequent analyses, 
we restrict our attention to the range 2 GeV^ < t < 4 GeV*. 

An analysis of the k = —1 normalized Gaussian sum-rule through ( |58| ) leads to predictions that are completely 
unstable under QCD uncertainties. Incorporating the error bounds of the dimension four giuonic condensate (52) and 
allowing for a 15% error in each of Hc and pc leads to a huge degree of variation in the resulting hadronic parameter 
estimates, and it is not even possible to ascertain whether or not an extension of the single resonance analysis is 
warranted. Mass predictions range anywhere from 1.0 GeV to 1.8 GeV, an interval far too broad to be of much use. 
The only inference we can draw from k — —1 analysis is that the mass scale obtained is in rough agreement with 
that which results from an analysis of the next-to- leading order (fc — 0) Gaussian sum- rule (see below). We note 
further that this consistency in mass scales between the k = —1 and higher order sum-rules is also observed in ||^, ^ 
and occurs only when instanton effects are included. Due to its extreme sensitivity to small variations in various 
QCD parameters, we are forced to conclude that the k — ~1 Gaussian sum-rule is an unreliable probe of the scalar 
glueball sector, and so, move on to a A; = analysis. 

A single narrow resonance analysis of the fc = Gaussian sum-rule, however, leads to rather poor agreement 
between theory and phenomenology. Minimization of ( |62|) leads to an optimum threshold at sq — 2.3 GeV^ which, 
when substituted into (63), yields a resonance mass of m = 1.30 GeV. In Figure ||, we plot both (Tq(t, sq) and 2r 
versus r for sq — 2.3 GeV^. The graph of (Jq{t, sq) appears to have a slope of two, but, if extended to r = 0, would 
not pass through the origin — a situation which contradicts the single resonance result (|65|). Further, in Figure ^, we 
plot the left- (theoretical) and right- (phenomenological) hand sides of ( psf ) versus s for r e {2,3,4} GeV* using the 
optimized values m = 1.30 GeV and Sq = 2.3 GeV^. The discrepancy between theory and phenomenology is apparent: 
the theoretical curves consistently underestimate phenomenology near the peak and overestimate phenomenology in 
the tails. These observations indicate that a single narrow resonance ansatz is an inadequate description of the scalar 
gluonium spectral function, hence motivating our subsequent analyses of more general models. 



3.2 Distributed Resonance Strength Models 

The moment combination measures the width of the QCD distribution, and since aQ — 2t > 0, we conclude that 
the resonance strength is distributed over an energy region broad enough to be sampled by the resolution •\/2t of 
the Gaussian kernel. In accordance, we should therefore consider phenomenological models which allow for such a 
spreading. We focus on three distinct classes of distributed strength models: 

1. single non-zero width models, 

2. a double narrow resonance model, and 

3. single narrow resonance + single non-zero width resonance models, 

ordered according to increasing complexity, i.e. the (normalized) single non-zero width models each contain two free 
parameters whereas the single narrow resonance -I- single non-zero width models requires four. In the subsections to 
follow, we outline the analysis procedure and the resulting hadronic parameter predications corresponding to each 
of the models under consideration. 



3.2.1 Single Non-Zero Width Models 

The first non-zero width model we consider is a unit-area square pulse representing a broad, structureless background 
which has previously been used for the study of the a resonance . To describe such a square pulse centred at 
and having width 2mr, we define 

-p^^<i(t) = -p'^P(t) = -^\e(t-m^+mr)~e(t-m^-mr)] , (66) 

TT TT 2ml 
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which leads to the foUowing normahzed Gaussian sum-rule (bTl) for k = 



1 



eri ■= — err 



2V7 



(67) 



As in the single narrow resonance analysis, the phenomcnological side of (|6^) admits a single r-independent peak, 
and so we can determine an optimum threshold parameter sq by following the steps outlined through (|60|)-(|6^) i.e. 
So = 2.3 GeV^. We wish to use this optimized value of sq to generate predictions for the hadronic parameters {m, F}. 
The following combinations of the moments 
in this regard: 



M) (where we suppress the explicit dependence on t and sq) are useful 



fe,i 



M, 



k.2 



A 



(3) _ Mk,3 „ fMk,, 



Mk,o \Mkfl) \Mk,o 



Mk,o 

MkA 



M. 



k,l 



(68) 



In the spirit of the sum-rules approach, we equate each of these theoretical quantities to its corresponding phe- 
nomcnological counterpart resulting in the following system of equations: 



Mi 



0,1 



M. 



0,0 



A, 



(3) 







(69) 

(70) 
(71) 



Considering that the right-hand (phenomcnological) sides of (|69|)-(|71|) are all r-independent, it is an important feature 
of the analysis that, at the optimized continuum threshold sq, the left-hand (theoretical) sides exhibit negligible 
dependence on r and so are well approximated throughout the interval [ti, Tf] by averaged values. Such r-independence 
of appropriate residual moment combinations also occurs in the more complicated models considered below. 

Then, we can use (|69|)-(|70[) to obtain predictions for the two free parameters {m, T} of the square pulse model (|6€). 
Equation ( [7l| ) is independent of our parameter estimates and therefore serves as an a posteriori consistency test!^ 

Inverting (||)-(^, we obtain m = (1.30±0.17) GeV and T = (0.59±0.07) GeV, where, in this and in all subse- 
quent analyses, the uncertainties quoted stem from the error bounds on the dimension- four gluonic condensate (^) 
and an estimated 15% uncertainty in each of ric and pc (see (|56|)). In Figure |[ we use sq — 2.3 GeV^ and central 
values of m and F to plot the theoretical and phenomcnological sides of ( [67| ) versus s for r g {2,3,4} GeV*. The 
excellent agreement between theory and phenomenology demonstrated by these plots is a vast improvement over the 
results of the k = single narrow resonance analysis (see Figure ^) . 

However, closer quantitative scrutiny reveals that the fits depicted in Figure are not quite as accurate as one 
might think. The quantity involves higher order moments of the fc = Gaussian sum-rule and so is sensitive to an 

even finer level of detail than are Aio,i/Mo,o and ct^. A QCD calculation of this quantity yields A'"q^ = —0.0825 GeV^ 
in significant violation of our consistency check (|7l|). 

For our second example of a single non-zero width model, we would ideally like to consider a Breit-Wigner 
resonance. However, closed-form expressions do not exist for the Gaussian image of a Breit-Wigner shape, and so 
we instead employ a Gaussian model to describe a well-defined resonance peak with a non-zero width 



9'^^^(i)=pS(t)^/2exp 



2F2 



(72) 



®For the square pulse model, it, se ems we could use either (^) or ( ^j] ) to compute the mass m; however, as we shall see, both yield the 
same result. We concentrate on ( |69[ ) simply because the analysis techniques outlined here are easily generalized to the more complicated 
resonance models to follow. 
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The quantity is a normalization constant related to the total resonance strength, and the Gaussian width T is 
related to an equivalent Breit-Wigner width Tbw through 

■ r 



Tbw = V21og2- 



(73) 



In this model, the k = normalized Gaussian sum-rule fel) becomes 



^o'^°°(s,T,so) 



where 



1 + erf 








V27r\/r2 + 


2t 







exp 



2 (r2 + 2t) 



(74) 



erf {x) 



dy 



(75) 



As in prior analyses, we begin with an extraction of the optimized continuum threshold parameter sq, and 
again, we accomplish this by examining the behaviour or the peak position. The phenomenological side of ( fz^ ) 
admits a single peak; however, in contrast to our previous two analyses, the position of this peak is dependent on r. 
Differentiating the right-hand side of (^3) and setting the result to zero gives 



TO 



exp 



' 4r^T(r^-|-2r) 



^ + °"V2rv^vT5+27 



(76) 



which, unfortunately cannot be explicitly solved for s. Consequently, in the absence of an exact solution, we approx- 
imate the phenomenological peak position by the expression 



A^ V — 



(77) 



where {A, i3,C} are to be considered unknown parameters. Explicit numerical experiments in (realistic) worst-case 
scenarios show that, provided r > 2 GeV*, the next term in the expansion (|77| ) [i.e. D/t'^] is negligible and can safely 
be ignored. Therefore, we are led to define the following x^.f^nction as a measure of the deviation of the theoretical 
peak position (Eoh from the phenomenological peak position characterized by the expansion (|77|): 



N r 



n=0 



Spoak(T-„, So) - A 



B 



C_ 

^2 



(78) 



The minimizing conditions 
dA dB dC 

can then be used to write {A, _B, C} as functions^ of sq leaving us with a one-dimensional minimization problem in 
So: 



(79) 



N 



Spcak(T„, So) - ^(So) - 



B{so) C(so) 



(80) 



'^The set of equations defined by condition ([i"9|) is linear and inhomogeneous. While trivial to obtain, the solution is rather a mess 
and so is omitted for brevity. 
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Minimizing ( |80| ) with respect to sq furnishes us with an optimized choice for the continuum threshold parameter. 

Once the optimized sq has been determined, the subsequent analysis of the Gaussian resonance model proceeds 
in a fashion completely analogous to the analysis of the square pulse model. For instance, the appropriate Gaussian 
model equations corresponding to (|69|)-(l7l|) are 



Mo,i 
Mo.o 



TA 



- 2r = - m^TA - T^A^ 



A, 



(3) 



r^) A + Sm^r^A^ + 2r^A 



^3 A 3 



where 



A 



exp 



(81) 

(82) 
(83) 

(84) 



Again, the first two equations may be (numerically) inverted to yield predictions for the model's two free parameters 
{to, r} while the third serves as an independent consistency check. The quantity A is small, and hence the terms 
in (^-|83|) proportional to A have a negligible effect on the extracted resonance parameters, a property which 
also simplifies our subsequent analysis. Note that since Aq'^^ ^ A in (|3|), we anticipate that this model will still 
underestimate the QCD value of this moment combination, a result which is confirmed in our detailed analysis. 

Carrying out the appropriate sequence of calculations first gives an optimized continuum threshold parameter^ 
So — 2.3 GeV^ which, in turn, yields a mass of to = (1.30 ± 0.17) GeV and an equivalent Breit-Wigner width 
(see ([73|)) of Fbvv = (0.40 ± 0.05) GeV. Again, we test the validity of these values by plotting both the theoretical 
and phcnomcnological sides of ( |7^ ) versus s for r S {2,3,4} GeV^, setting sq — 2.3 GeV^ and employing the 
central values of the resonance parameters. The resulting graphs are, to the naked eye, indistinguishable from those 
corresponding to the square pulse model analysis^ (Figure |) and so are omitted. 

Substitution of the central values of to and F into (^ ) yields A^^"^ = 0.000342 GeV^ which must again be compared 

with the QCD value of ^q^' ~ —0.0825 GeV^. Clearly, there exists a significant discrepancy between the two as the 
Gaussian model cannot even correctly predict the sign of this moment combination. 

The quantity Aq^-*, however, is a measure of the asymmetry of N^'^^{s, r, sq) with respect to s about its average 
value defined by Mo,i/Mo,o- Therefore, considering that both the Gaussian and square pulse models represent reso- 
nance strength distributions that are symmetric (about to^), it is perhaps not surprising that they fail to accurately 



(31 

predict Aq . Hence, we are prompted to consider a skewed generalization of the Gaussian resonance model: 



{t - m^ f 
2f2" 



(85) 



where the factor introduces a degree of asymmetry and has been chosen to achieve consistency with known 
low-energy two pion decay rates [|lO[ ^ . Substituting (|5|) into ( |5l| ) gives 

(s — m?Y "' 



£;A^o (s.^>so) = 2V¥exp 



+ 4exp 



2(F2 + 2r) 



{fV^ + AsT'-m'-T + AiJi^T^ + 2rF'* + 4t^F^) 



4_2 



1 + erf 



rV^x/r^ + 2t(sf2 + 



where 



2TO^Fexp(-|J') V^(r' + 2rf/' 



V2FV7VF2 + 2t 
(86) 
(87) 



Surprisingly, this is the same optimized threshold narametor value that we found in the single narrow resonance model analysis. This 
need not always be the case as depicted explicitly in Bl]^ 
®We quantitatively address this situation in Section W 
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The sum-rules analysis proceeds along the same lines as for the unskewed Gaussian resonance model. The appropriate 
generalizations of (pl|)-(p3|) are 



Mo,i _ m2(m4 + 3r2) 



a: 



2r 



m4 + r2 

(m4 + r2)2 
(3) _ Am^r^{m^ - 3r2) 

(TO4+r2)3 



0(A) 
+ 0(A) 
+ 0(A) 



(88) 
(89) 
(90) 



Numerically inverting[J (^8|)-(|89|), we find a mass m = (1.17 ± 0.15) GeV and an equivalent Breit-Wigner width of 
Tbw — (0.49 ±0.06) GeV. A subsequent phenomenological prediction yields ^q^'' — 0.00943 GeV^ — again, completely 
inaccurate compared with the QCD value. These results, as well as the results of the unskewed non-zero resonance 
width model analyses are summarized for convenience in Table 1. 



resonance model 


mass (GeV) 


width (GeV) 


(GeV6) 


unskewed Gaussian 


1.30 


0.38 


0.000342 


unskewed square pulse 


1.30 


0.59 





skewed Gaussian 


1.18 


0.49 


0.00943 



Table 1: The results of a /c = Gaussian sum-rules analysis of a variety of non-zero resonance width models using 
central values of the QCD parameters. For the Gaussian resonance models, the given width is actually the equivalent 



Breit-Wigner width (see ([73|)). The A^^"^ values should be compared with the QCD prediction A, 



(3) 



-0.0825 GeV^ 



3.2.2 Double Narrow Resonance Model 

The double narrow resonance model is defined by 



(91) 



where mi < m2 are the two resonance masses and /i,/2 are their respective couplings. Correspondingly, the 
normalized Gaussian sum- rule (|l]) reduces to 



exp 



(s — rrilY 



f^mf 



: exp 



At 



(92) 



In the analysis of this model, it becomes inconvenient to use the parameters {toi, /i, ?n2, /2} and so we instead focus 
on the set {z, y, r} defined by 



n - r2 



with 



/2mf 



ri 



r2 = 



with ri + r2 = 1 



(93) 



(94) 



As in the Gaussian resonance model analysis, the phenomenological side of (|9^) admits a single peak whose 
position is r-dependent. In terms of the double resonance model parameters (p3|), differentiating the right-hand side 
of (p2) with respect to s and setting the result to zero yields 



(^„l)(s-iz+iy) 



exp 



viz -2s) 



At 



= 



(95) 



"The terms in 



891) proportional to A arc found to have a negligible impact on the solutions. 
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which, again, cannot be explicitly solved for s. Thus, we approximate the phenomenological peak position by the 
expression ( [tT] ) and correspondingly, extract an optimized continuum threshold parameter sq = 2.3 GeV^. 

To compute predictions for the hadronic parameters of the double narrow resonances model, we again look to 
moment combinations such as (^8|). However, since this model contains three free parameters (p3[), we require three 
equations. In other words, if we wish to have an independent consistency check of our results, we must introduce a 
fourth moment combination in addition to those of (p8|). Defining 

(96) 



Sk = 



Ml 



kA 



,Mk,3Mk,i , ,Mk,2 fMk,i 



Mkfl Mk,o Mk,o 
leads to the following equations 
Mfci 1 



6- 



M, 



k,0 



fe.O 



-3 



fe,i 



fc,0 



M, 



k,a 



-(z + ry) 



2r=-y'(l- 



l(3) 



rV(i 



1 



12T{al--2T)^-y^l-r') {l + 3r') 



(97) 

(98) 
(99) 

(100) 
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which are analog ous to (ii)-(0) and tel^)-(l83). 

For the fc = case of interest, upon inverting (|97|)-(p9|) we find the heavier of the two states is also the more 
strongly coupled with m2 = (1.4 ± 0.2) GeV and r2 = 0.72 ± 0.06. The lighter resonance with a mass mi = 
(0.98 ± 0.2) GeV is the more weakly coupled state with ri = 0.28 =F 0.06. The QCD uncertainties given for the 
resonance masses obscure the very stable mass splitting between the two states: m2 — nii = (0.42 ± 0.03) GeV. 
Again, plots of the theoretical and phenomenological Gaussian sum-rules exhibit the excellent agreement shown in 
Figure ||. 

Using central values of the hadronic parameters obtained, we can calculate a phenomenological prediction associ- 
ated with the first independent moment combination (|9^). The result is Sq - 12t^ - 12t {a^ - 2t) = 0.073733 GeV^ 
which must compared with a QCD calculation of - Ut'^ - 12t (ag - 2r) = 0.169769 GeV^. The roughly 50% de- 
viation between these values is a large improvement over similar comparisons in the single non-zero width resonance 
models. 



3.2.3 Two Resonance Model of a Narrow Resonance Combined with a Wide Resonance 



A natural extension of the double narrow resonance model is to introduce an additional parameter into the hadronic 
model which describes the width for one of the resonances. 

We begin with a model consisting of a narrow resonance of mass m and a Gaussian resonance of mass M and 
width r. The resulting k = normalized Gaussian sum- rule (|l|) is: 



1 



; exp 



(s — 



At 



tm exp 



'2(r2-H2r) 



1 + erf f sr^+2M^r 



27rVr2 + 2t 



1 -I- erf 



V2r 



where and vm denote the relative strengths of the resonances and are constrained by r„ 



1. 



(101) 
As in the 



other models, the phenomenological side of (101) has a single r-dependent peak position. The peak-drift x" in (|7q) 



is again minimized to find the optimum value of the continuum, and then the theoretical (QCD) values for the 
moments are used to determine the resonance parameters. Since this model contains four independent parameters, 
the lowest four moment combinations (|6^ , |96| ) are used to determine the parameters and the following fifth order 
moment combination will be used as a consistency check. 

2 ,^ /,^\3 /,^-\5 



l(5) 



Mk' 



.MkAMk; 



Mk,o Mkfi Mkfi 



10 



Mk,3 f Mkj 



Mk.o \Mk.o 



10 



Mk^2 f MkA 



Mkfi \Mkn 



Mk.o 



(102) 
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Defining 

z = + m? , y = rr? — Af^ , r = r„j — rjv/ 
we find the following expressions for the moment combinations in terms of the resonance parameters 
Mo,i 



1 



0,0 



(z + ry) + O (A) 



ag-2r=i2;2(l-r2) + V(l-r)+0(A) 



l(3) _ 



0(A) 



S, - I2r^ - 12r (a^ - 2r) = ly^ „ ^2) + 3^2) ^ 3^2^2 ^ ^) _ ^2) 



(103) 

(104) 

(105) 
(106) 
(107) 
(108) 



The terms proportional to A (see equation (p^)) in the above expressions are found to be numerically insignificant. 

The QCD values of the moments are used to solve (104- 107[ ) for the resonance parameters.]^ The non- linear 
nature of (104-107) might suggest a large number of solutions, but it is found that only one physical solution 
persists across the range of QCD parameters considered: m — (1.41 ± 0.19) GeV, M — (1.23 ± 0.15) GeV, Tbw = 



(0.52 ± 0.06) GeV, r,„ = 0.49 ± 0.13, tm = 1 



Plots of the theoretical and phenomenological sides of the 



normalized Gaussian sum-rules are again represented in Figure |^, illustrating excellent agreement between the QCD 
prediction and phenomenological model. 

The fifth-order moment (108) can now be used as a consistency check. The resonance parameters for the central 



QCD parameters yields A^^^ - 20r^*,^^ 
20tA[,^^ = -0.243GeVi'^. 



0.130 GeV which should be compared with the QCD value of A, 



To assess whether other non-zero width models lead to better agreement with the QCD value of Aq^-* — I^StJ^^^ = 
— 0.243 GeV^° we consider using the squ are pulse model instead of the Gaussian resonance, which results in a 
normalized Gaussian sum- rule similar to ( 101 ) 



(5) 



QCD/- N 

\s,t,bq) = t^ 



'Attt 



: exp 



n2\2 



4t 



1 



rM 



AMY 



erf 



s- AP + MT 



erf 



AP - MT 
2^ 



(109) 



where -f ~ 1. Using the definitions (103) results in the following expressions for the moment combinations in 
terms of the resonance parameters. 



Mo,: 
M, 



0,0 



(3) 



i (z + ry) + O (A) 
2r = _ ^2) ^ _^p2 _ (1 _ ^) + o (A) 

-\v\ (1 - r^) - Ir^y (1 - r^) (z - y) + O (A) 



5o - 12t2 - 12t {gI - 2r) =^^y^ (l - r^) (l + Zr^) + ip^y^ (1 + r) (l - r^) (z - y) 



a; 



(5) 



20tA 



(3) 



+ lr4(l-r) (z-y)VO(A) 
\y\ (1 - r^) (1 + r^) - ^F^y^ (l - r^) (1 + rf (z - y) 
Ir4y(l-r2) (z-y)VO(A) 



(110) 

(111) 
(112) 

(113) 
(114) 



^^As mentioned earlier, the OC p v a lue of the residual moment combinations are found to be virtually t independent consistent with 
the resonance model expressions (L04-L03). 
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Using the QCD values of the moments to solve ( |110 -113) for the resonance parameters again yields a single 
physical solution over the QCD parameter space: m = (1.33 ± 0.18) GeV, M = (1.23 ± 0.18) GeV, Tbw = (0.95 ± 
0.12) GeV, Tm = 0.60 ± 0.13, tm = 1 — r^- Agreement between the theoretical and phenomenological sides of the 
normalized Gaussian sum- rules is again excellent, as exhibited by the plots shown in Figure ||. 

The value of the fifth-order moment combination (114) for the resonance parameters corresponding to the central 
QCD parameters is A^^ — 20tAq'^' ~ — 0.114GeV^°. The square pulse model is thus not an improvement upon the 
Gaussian resonance model's agreement with the QCD value of this resid ual m oment combination. 

As a final scenario, we modify the normalized Gaussian sum-rule of (101) for a skewed Gaussian resonance 




(s — -mr) 



2\2-\ 



At 
47f2 



1 + erf 



/ sr2 + 2rA/2 



rV7\/r2 + 2t{sT'^ + 2tM^) 



\2Ty/^^/r^ + 2t 
(115) 



where E (with appropriate substitution of the mass M) is defined in (p7|). Using the same conventions, we obtain 
lengthy expressions for the moment combinations in terms of the resonance parameters. 



D 



Mo- 



Mo,o 
D = {z- 



1 



(z + ry) {z - yf + {-2y + 4z + Ary - 2rz) + O (A) 



y)^ + 4r^ 



(116) 
(117) 



{al - 2r) =-v^ (l - r^) [z - yf + -V\l - r) {z - yf {z^ - Aryz - 6yz + + gry^) ^^^^^ 
+ AT^ (1 - r^) [z - 2y f + 2AT^{1 - r) + O (A) 

D^AI^^ = - iyV (1 - r^) (z - yf - ^T^y (l - r^) (z - yf (z^ - 4zry - 2yz + 8ry^ + y^) 
+ (l - r^) (z - yf (z - 2y) (z^ - Azry - 2yz + 8ry'^ + y^) 

+ AT'^il - r) [y^ (-32r^ - 41r - 13) + z-''4 (r^ + r + l) + y^z6 (8r^ + llr + 5) + z^y {-2Ar^ - 33r - 2l)] 
+ 48r^(l - r) {~6ry + 'irz - 2y - z) + O (A) 

(119) 
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[5o - 12r2 - 12r {al 2r)] = (l - r^) (l + Sr^) (z - y)« 

+ ir^y^ (l - r^) (z - yf {irz^ + 3z^ - Qzry - Uzr'^y - IQyz + 2Ar^y^ 



lly^ + W) 



+ -r* (1 -r){z- yf {z^ + y"- [A&r' + SGr^ + 32r + 25) + (-4r2 - 8r - 8) + j/^z (-48r^ - GSr^ 

(I2r3 + 28r2 + 36r + 26)] 



56r 



+ (1 - r^) (z - y)^ [z43(l + r) + y4 (gOr^ + 21r + 6l) + yz^ {-Vlr^ - 18r - 30) + y^z (-144r2 - 54r - 1 

(72r2 + 48r + 96)] 

+ 16r^ (1 - r) [z^ (3r^ + 'ir^ + 9r + 12) + (48r''' + 84r2 + 80r + 47) + y^z (-OGr^ - ISGr^ - 196r - 118) 

+z3y (-24r^ - 42r2 - 60r - 54) + y^z^ {72r^ + 144r^ + 168r + 114)] 

+ 192ri° (1 - r) [y^ {Ur^ + 8r + 8) + z^ (3r2 - 2r + 7) - 12yz(l + r)] 
+ 192Ori2(-^_j.)^0(^) 

(120) 



(5) 



20rA 



(3) 



Vr(l-r2)(l + r2)(z-2/) 
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- -r^ (1 - r^) (z - yf z^ (1 + r)^ + y^ (Sr^ + + lOr + l) + yz {-4r^ - 2r^ - 8r - 2) 
8 L 

- ^r*y (l - r^) (z - y)^ [3z'' + y* (64r^ + 12r2 + 88r + 15) + z^y {-6r^ - 12r - 18) 

+y^z (-64r3 - 30r2 - 124r - 42) + y^z^ {l6r^ + 2Ar'^ + 64r + 42)] 

+ bV^ (l - r^) (z - yf [3z^ + y^ (-128r^ - 30r2 - 196r - 44) + y^z (l92r^ + 84r2 + 384r + 135) 
+ z^y (-6r2 - 12r - 24) + z^y"^ (iGr^ + 36r2 + 96r + 86) 
+y^z^ (-96r3 - 84r2 - 288r - 156)] 

+ 40r^ (1 - r) (z - yf [y" (-64r'' - 90r^ - 146r2 - 153r - 37) + z'^ (r^ + ir"^ + 3r + 5) 



36r2 - 47r - 39) 



+ y^z {l2Sr^ + 201r3 + 367r2 + 399r + 125) + z^y {-Ar^ - 12r3 
+ z^y^ (-327,4 ^ ggj,3 _^ -^g2r2 + 202r + 114) 

+y^z2 (-ger-^ - 170r3 - 354r2 - 408r - 168)] 
32ri°(l - r) [y^ (-128r'* - 308r3 - 508r2 - 558r - 194) + z^ {Ar'^ + Ar"^ + 24r2 + 59r - l) 

+ z^y {-AQr^ - 85r3 - 195r2 - 485r - 155) + y^z (320r'* + 860r3 + 1420r2 + 1755r + 695) 
+z3y2 (I60r'' + 430r3 + 750r2 + 1460r + 620) + y^z^ (-320r'' - 905r3 - 1495r2 - 2235r - 

.3 tn^Z 



+ 960r^2(l - r) [z^ {2r^ - 2r^ + 14r - 6) 



-16r^ - 16r^ - 37r - 13) 



+z2y (-12r3 + Ar^ - 57r - l) + y^z (24?'^ + 8r2 + 74r + 18)] 



+ 3840r"(l - r) [y(-10r - 2) + z(5r - 3)] + O (A) 



(121) 



Using the QCD values of the moments to solve (116-12C) for the resonance parameters again yields a single, 
stable, physical solution over the QCD parameter space: m = (1.38 ± 0.13) GeV, M = (1.06 ± 0.21) GeV, T bw = 
(0.69±0.07) GeV, r™ = 0.44±0.04, tm = 1 — ?'m- Excellent agreement between the theoretical and phenomenological 
sides of the normalized Gaussian sum-rules is illustrated by the plots shown in Figure ^ 

The value of the fifth-order moment (121) for the resonance parameters corresponding to the central QCD 



parameters is A\^'' 



20rA 



(3) 



-0.192 GeV^", a result which has only a 21% deviation from the QCD value AJ,' 



(5) 
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Resonance Model 


m (GeV) 


M (GeV) 


Width (GeV) 


r = r.m - rM 


Double Narrow 


1.41 


1.00 





0.42 


Narrow plus Square 


1.34 


1.24 


0.93 


0.18 


Narrow plus Gaussian 


1.38 


1.23 


0.51 


-0.10 


Narrow plus Skewed Gaussian 


1.40 


1.00 


0.68 


-0.16 



Table 2: Resonance parameters obtained from central values of the QCD parameters in the various two-resonance 
scenarios. The mass M denotes the state associated with the quoted width. The width parameter quoted for the 
Gaussian models is the equivalent Breit-Wigner width. 





So - 12t2 - 12r (ct2 - 2t) (GeV^) 


4'^-20rA^^^ (GcVi") 


QCD 


0.1698 


-0.2433 


Double Narrow 
Narrow plus Square 
Narrow plus Gaussian 
Narrow plus Skewed Gaussian 


0.0737 


-0.0502 
-0.1144 
-0.1300 
-0.1918 



Table 3: Comparison of the next-highest moment combinations in the two-resonance scenarios with the QCD values. 
All values are obtained from the central values of the QCD parameters. 

20tAI^'' — — 0.243 GeV^". The skewed Gaussian model is thus in reasonable agreement with the QCD value of the 
fifth-order residual moment combination. 

4 Results and Discussion 

The resonance parameters devolving from the central values of the QCD parameters have already been summarized 
m Table I for the single wide resonance models, and Table H summarizes the results for the resonance parameters 
obtained in the double resonance models. We note that in the narrow plus wide resonance scenarios, no assumption 
was made on which state would be wide, and hence it is interesting that the analysis consistently predicts that the 
lightest state has the non-zero width. Similarly, Tables |l| and ^ summarize the QCD and predicted values of the 
moments which serve as a consistency check in each of the models. Based on this criteria, the double resonance 
models have much better agreement with the QCD prediction than the single wide resonance models. The narrow 
plus wide resonance model is the most accurate in this respect, with only a 21% deviation from the QCD value of 
the fifth order asymmetry moment Aq^'' — 20tA'^^ . 

As illustrated in Figure ^ the single narrow resonance model leads to a significant deviation between the theo- 
retical (QCD) and phenomenological results for the normalized Gaussian sum-rule. The single wide resonance and 
double resonance models dramatically improve this agreement, and all the models lead to plots which are represented 
by Figure |5| However, a measuring the difference between the theoretical and phenomenological curves over the 
ranges of s and r used in Figure |^ can be used as another criteria to evaluate the effectiveness of the phenomenological 
models since it provides a quantitative measure of the difference between the theoretically and phenomenologically 
determined fc = Gaussian sum-rules. Table ^ lists this the optimized value so — 2.3 GeV^ in each of the 

resonance models considered for the central QCD parameters. We see that the double resonance models provide a fit 
which is an order of magnitude better than any of the non-zero width resonance models. However, the is on the 
order of 10^® for all the double narrow resonance models, and hence does not provide a strong distinction between 
the various scenarios. 

Thus the moment consistency test and measure of the agreement between the phenomenological and QCD 
values of the Gaussian sum-rules clearly favour the double resonance scenarios. However, this does not clearly 
distinguish between the various two-resonance scenarios. Despite this apparent difficulty in distinguishing between 
the scenarios, the mass of the two states in the various double resonance scenarios are relatively stable, lying in 
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resonance model 




unskewed Gaussian 


1.21 X 10"'' 


unskewed square pulse 


1.31 X lO-'' 


skewed Gaussian 


1.45 X 10-5 


double narrow 


1.83 X 10"" 


narrow and square pulse 


0.923 X 10-t> 


narrow and unskewed Gaussian 


0.888 X 10-« 


narrow and skewed Gaussian 


0.959 X 10-« 



Table 4: The x for the fits between the theoretical and phcnomenological results for the normalized A; = Gaussian 
sum-rules over the ranges of s and r used in Figure || for each of the phcnomenological models considered. Central 
values of the QCD parameters have been employed. 

the range 1.0-1.4 GeV, with the greatest mass splittings occurring in the double narrow resonance model and in the 
narrow plus skewed Gaussian resonance model. 

5 Conclusions 

In this paper, we used QCD Gaussian sum-rules to analyze the scalar glueball sector in an effort to obtain predic- 
tions for the hadronic parameters {i.e. mass, width, coupling strength) of low-lying scalar glueball states. In our 
analysis, we incorporated instanton effects and employed a number of phcnomenological models more general than 
the traditional single narrow resonance. 

First, we demonstrated that the leading order k = —\ Gaussian sum-rule (the only sum- rule which receives a 
contribution from the LET) leads to results which are unstable under moderate QCD uncertainties: single resonance 
mass extractions ranged anywhere from 1.0 GeV to 1.8 GeV. This variation is too large to yield a definitive prediction 
for the mass of the lightest scalar glueball, but does, however, indicate a lower bound of roughly 1 GeV on gluonium 
mass scales. This observation supports similar results obtained previously in ||^, || where it was shown that instanton 
effects serve to increase the scale of masses extracted from the k = —1 sum-rule whereby reducing the discrepancy 
between masses predicted using this k = — \ sum-rule and those extracted from fc > sum-rules. 

Due to the instability of the fc = — 1 Gaussian sum-rule analysis, we turned to the fc = sum-rule. We showed 
that hadronic parameters extracted using a single narrow resonance model led to poor agreement between theory and 
phenomenology as indicated by the plots of Figures |^ and ^. Consequently, we then considered phcnomenological 
models which allowed for resonance strength to be distributed over an appreciable energy range. 

We focused on three single non-zero width resonance models (see Section ||), and the results of all the corresponding 
analyses were indicative of a (m « 1.2-1.3 GeV), wide (F k, 0.4-0.6 GeV) resonance (see Table The coincidence 
between plots of the theoretical and phcnomenological Gaussian sum-rules (characterized by Figure ^ represented 
a vast improvement over the corresponding graphs obtained from a single narrow resonance analysis (Figure U). 
However, as indicated in Table [U all three models failed completely to predict the value of the first independent 

(3) . 

moment combination , prompting us to move on to analyses of phcnomenological models that allow for a second 
resonance. 

Regarding both a measure of the discrepancy between theoretical and phcnomenological fc = Gaussian 
sum-rules and the ability to accurately predict higher order, independent moments, the double resonance models 
represented a marked improvement compared to the single non-zero width resonance models. Resulting hadronic 
parameter estimates from the various models analyzed are summarized in Table ^. Estimates of the larger of the 
two masses (m) were remarkably insensitive to the particular width model used, approximately 1.4 GeV in all cases 
considered. Mass predictions for the lighter of the two resonances (M) ranged from 1.0 GeV in the double narrow 
and narrow plus skewed Gaussian models to roughly 1.25 GeV in the narrow plus square and narrow plus Gaussian 
models — again, a surprising degree of stability considering the variety of width models employed. Further, in those 
models which allowed for a non-zero resonance width, we made no a priori assumptions as to which of the two 
resonances was wide; however, the width was consistently found to be associated with the lighter of the two states 
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and was relatively large in magnitude, yielding an equivalent Breit-Wigner width of Few ~ 0-55 GeV in the Gaussian 
resonance models. These hadronic parameter estimates were subsequently employed in predictions of higher order, 
independent moment combinations. The results of this analysis are summarized in Table ^. The deviation of the 
phenomenological predictions from the values computed using QCD was typically 20-50%, with the best agreement 
occurring in the skewed Gaussian resonance model. While certainly not perfect, this represented a vast improvement 
over analogous calculations performed in the various single non-zero width resonance models (see Table |^) . Finally, 
in Table ^, we have collected the results of a measure of the difference between the theoretical k = Gaussian 
sum-rules and the various phenomenological sum-rules corresponding to each of the resonance models considered. For 
the double resonance models, the values quoted were consistently an order of magnitude lower than those pertaining 
to the single non-zero width models. Therefore, in addition to our previous conclusion concerning the necessity 
for distributed resonance strength, we note further that a single non-zero width resonance model is insufficient to 
account for this spreading and that far better consistency between theory and phenomenology (as encapsulated in 
the Gaussian sum-rules of scalar gluonium) is achieved by employing models which allow for a second resonance. 

Although certainly relevant to the nature of the spectrum of scalar states with masses under 2 GeV, it is difficult 
to make a direct comparison between our results and the entries of the PDG since the typical mass separation 
between the various observed scalar states is relatively small: on the order of 0.2 GeV. Unfortunately, this is 
approximately the same magnitude as the uncertainties we calculate in our mass estimates. These error bars are the 
direct result of uncertainties in the QCD parameters ( 
from drawing definitive conclusions concerning the nature of any specific scalar state. Instead, we merely note that 
certain classification schemes proposed by Minkowski and Ochs and Narison Q require a relatively light (« 1 
GeV), wide state — a situation completely consistent with our results. Furthermore, the persistent m = 1.4 state we 
find may indeed have ramifications regarding the nature of either the /o(1370) or the /o(1500) p8[ |. 

Acknowledgements: The authors are grateful for research support from the Natural Sciences and Engineering 
Research Council of Canada (NSERC). 



Therefore, in light of these uncertainties, we refrain 
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Figure 1: Contour of integration Fi + F2 defining the Gaussian sum-rule in (Pq). The wavy hne on the negative real 
axis denotes the branch cut of n(w). 
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Figure 2: Closed contour C{R) used to calculate the Gaussian sum-rule defined by (p5|). The inner circular segment 
has a radius of e, and the circular segments F3, F4 and F5 have a radius R. The wavy line on the negative real 
axis denotes the branch cut of ll{w), and the linear segments of the contour above and below the branch cut are 
denoted by F^,. The contour f i(-R) is that portion of Fi (see Figure P which lies in the interior of a circle of radius 
R centered at — s, and similarly for F2(i?). 
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Figure 3: Plot of CTq for the theoretical prediction (dotted curve) compared with ctq = 2t for the single-resonance 
model (solid curve) for the k = sum-rule using the x^-optimized value of the continuum so = 2.3 GeV^. 
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Figure 4: Comparison of the theoretical prediction for Nq (s, t, sq) with the single narrow resonance phe- 
nomenological model (^8|) using the x^-optimized values of the resonance mass (m — 1.30 GeV) and continuum 
(so = 2.3 GeV^). The r values used for the three pairs of curves, from top to bottom in the figure, are respectively 
T — 2.0 GeV*, T — 3.0 GeV'', and r = 4.0 GeV*. The phenomenological model is consistently larger than the 
theoretical prediction near the peak, but is consistently smaller than the theoretical prediction in the tails. Central 
values of the QCD parameters have been used. 
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Figure 5: Typical comparison of the theoretical prediction N^^^ (s, r, sq) with the phenomenological side of 
the normalized Gaussian sum-rule (|l|) in the single non-zero width resonance models and the double resonance 
models. Individual plots for each model are indistinguishable from those shown, and are hence not repeated as 
additional Figures. The x^-optimized value of the continuum (sq = 2.3 GeV'^) has been used to extract the resonance 
parameters as outlined in the text. The r values used for the three pairs of curves, from top to bottom in the figure, 
are respectively t = 2.0 GeV^, r = 3.0 GeV'', and r = 4.0 GeV^. Note the almost perfect overlap between the 
theoretical prediction and the phenomenological models. Central values of the QCD parameters have been used. 
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